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Inclusive Standard Model Higgs boson pair production and sub- 
sequent decay to same-sign dileptons via weak gauge W bosons 
at the CERN Large Hadron Collider has the capability to determine 
the Higgs boson self-coupling, A. The large top quark mass limit is 
found not to be a good approximation for the signal if one wishes to 
utilize differential distributions in the analysis. We find that it should 
be possible at the LHC with design luminosity to establish that the 
Standard Model Higgs boson has a non-zero self-coupling and that 
A/Asm can be restricted to a range of 0-3.7 at 95% confidence level 
if its mass is between 150 and 200 GeV. 

The CERN Large Hadron Collider (LHC) is widely re- 
garded as capable of directly observing the agent responsible 
for electroweak symmetry breaking and fermion mass gen- 
eration. This is generally believed to be a light Higgs bo- 
son with mass M H < 200 GeV [[[)]. Furthermore, the LHC 
promises complete coverage of Higgs decay scenarios [Q], in- 
cluding general parameterizations in the Minimal Supersym- 
metric Standard Model invisible Higgs decays |Q], and 
possibly even Higgs boson decays to muons [gp. This broad 
capability was made possible largely by the addition of the 
weak boson fusion production channel to the search strate- 
gies [|3|JJ. Observation of a Higgs boson in this channel also 
contains additional information in the angular distributions of 
the scattered quarks which reveal the fundamental tensor struc- 
ture of the VVH vertex [Q]. With mild theoretical assump- 
tions, the Higgs boson total width, Th, can be determined 
via combination of all available channels, which in turn yields 
the gauge and various Yukawa couplings [^J^j. The weak bo- 
son fusion channels have received considerable attention in the 
LHC experimental collaborations, and a number of detailed 
detector simulation studies on them have already been com- 
pleted, with very encouraging results Jl0[ ] . 

While these studies have shown that the LHC promises 
broad and significant capability to measure various properties 
of the Higgs sector, what remains is to determine the actual 
Higgs potential. This appears in the Lagrangian as 

V($) = -Au 2 ($t$) + A($t$) 2 , 

where $ is the Higgs field, v = (v^Gf) -1 ^ 2 is the vacuum 
expectation value, and Gf is the Fermi constant. In the Stan- 
dard Model (SM), A = X S M = Mfj/(2v 2 ). Regarding the 
SM as an effective theory, A is per se a free parameter. Its up- 
per limit can be determined using unitarity arguments, assum- 
ing the model's validity to high energy scales [11]. To measure 
A, and thus determine the Higgs potential, at a minimum ex- 



periments must observe Higgs boson pair production; while 
this has been shown to be possible for a light Higgs boson at 
a future Linear Collider [ 12], no study has yet been presented 
which suggests this is possible at the LHC for the SM Higgs 
boson. 

We simulate the signal process, pair production of two SM 
Higgs bosons in gluon fusion, at the parton level for pp colli- 
sions at ^/s = 14 TeV. Both Higgs bosons are decayed to W 
boson pairs, which subsequently are decayed to four jets and 
two same-sign leptons [][^| 
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HH -> (W + W~) (W + W~) -> (jj^v) (jjl v), 



where £, I' are any combination of electrons or muons. The 
intermediate Higgs and W bosons are treated off-shell using 
finite widths in the double pole approximation. We calculate 
the signal using two methods: exact loop matrix elements [ |l4| ] 
and the infinite top quark mass limit. The latter, which is com- 
monly used in place of exact matrix elements to speed up the 
calculation, reproduces the correct total cross section for HH 
production to within 10% to 30% for Higgs masses between 
140 GeV and 200 GeV. However, it produces completely in- 
correct kinematic distributions. 

Signal results are computed consistently to leading order 
QCD with the top quark mass set to m t — 175 GeV and 
SM top Yukawa coupling, and the renormalization and fac- 
torization scales are taken to be the Higgs boson mass [|4D. 
The effects of next-to-leading (NLO) order QCD corrections 
are included by multiplying the differential cross section by an 
overall factor K = 1.65 (if-factor), as sugg ested by Ref. 
where the QCD corrections for gg — > HH have been com- 
puted in the large m t limit. The multiplicative effects of NLO 
QCD corrections are not expected to depend on whether the 
signal is calculated with exact matrix elements or in the infi- 
nite top quark mass limit. 

The SM backgrounds of interest are those that produce two 
same-sign leptons and four well-separated jets which recon- 
struct in two pairs to a window around the W boson mass. The 



1 While our study was in progress, Ref. jl3| ] appeared. It includes a 
brief discussion of Higgs boson pair production at an upgraded LHC, 
which would gather 20 times the amount of data expected in the first 



run. 
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Unfortunately the search in this channel cannot be generalized to 
the supersymmetric case, since the branching fraction to W bosons is 
suppressed for the light Higgs scalar close to the decoupling regime. 
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Figure 1. Minimum separation between jets, AR(jj)min, for the 
Mh = 180 GeV, signal with exact matrix elements (solid line) and 
in the large rot limit (dot-dashed line), and the WWWjj (dashed 
line) and tiW backgrounds (dotted line). Qualitatively similar results 
are obtained for other values of the Higgs boson mass in the range 
150 GeV <M H < 200 GeV. 

largest contribution originates from W ± W + W~ jj produc- 
tion, followed by ttW ± where one top quark decays leptoni- 
cally, the other hadronically, and neither b quark jet is tagged. 
Other backgrounds, which in sum contribute at the < 5% 
level Jl3|], are: tttt production, where none of the b quark 
jets are tagged, and additional jets or leptons are not observed; 
W Zjjjj production with leptonic Z decay (including off- 
shell photon interference) where one lepton is not observed; 
and ttj events where one b quark decays semileptonically with 
good hadronic isolation, and the other is not tagged. For this 
letter we consider only the two dominant backgrounds; since 
the others enter in sum at less than 5% of the total contribution, 
they do not change our results noticeably. 

We simulate both leading backgrounds at the parton level 
using exact matrix elements generated with MADGRAPH Qlq]. 
For WWWjj production we evaluate the strong coupling 
constant a s and the parton distribution functions at a scale /i 
given by /.i 2 = J2Pt> wnere the sum extends over all final 
state particles; for ttW production we take fi = m t + My/ /2. 
We use a value for the strong coupling constant of a s (Mz) = 
0.1185. QCD corrections are not taken into account in our cal- 
culation of WWWjj and ttW production. The top quarks are 
generated on shell (narrow width approximation), while all W 
bosons in both processes are allowed to be off shell. Assuming 
a b quark tagging efficiency of 50%, only 1/4 of the ttW rate 
contributes to the background; events with one or two tagged b 
quarks are rejected. All signal and background cross sections 
are calculated using CTEQ4L [ |l7| ] parton distribution func- 
tions. 

The kinematic acceptance cuts for both signal and back- 
grounds are: 
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Figure 2. Distribution of the invariant mass of all observ- 
able final state particles, m vis , after all cuts, for the signal with a) 
Mh = 160 GeV and b) Mh = 200 GeV, and the dominant back- 
grounds. The m V i 3 distribution of the signal evaluated in the large 
ra t limit is also shown. 

PT (j) > 30, 30, 20, 20 GeV, p T (i) > 15, 10 GeV, 

IrKj)! < 3.0, W)\ < 2.5, 

AR(jj) > 0.6, AR(j£) > 0.4, AR(U) > 0.2, 



{A(j}) + (Any 



1/2 



is the separation in the 



where AR = 

pseudorapidity - azimuthal angle plane. In addition we re- 
quire the four jets to combine into two pseudo-T'K pairs with 
invariant masses between 50 and 110 GeV and assume that 
this captures 100% of the signal and backgrounds. We do not 
impose a missing transverse momentum cut which would re- 
move a considerable fraction of the signal events. 

Both backgrounds are multi body production processes, and 
one therefore expects that the distribution of the invariant 
mass, \fl, of the system peaks at values significantly above 
threshold. In contrast, the signal is a two-body production 
process for which the \f§ distribution will exhibit a sharper 
threshold behavior. Since the Higgs bosons are produced 
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Figure 3. The m vis distribution of the signal for Mh = 180 GeV 
in the SM (solid curve), for Xhhh = A/Asa/ = (dashed line) and 
for Xhhh = 2 (dotted line). The dot-dashed line shows the SM 
cross section in the large m t limit. Qualitatively similar results are 
obtained for other values of Mh ■ 

almost at rest, the final state particles are distributed fairly 
isotropically, resulting in a distribution of the minimum jet-jet 
separation, Ai?(jj) min , which is peaked at AR(jj) min w 1.3 
(see Fig. [j]). In contrast, the Ai?(jj) m ; n distribution for the 
background processes peaks at at a lower value (AR(jj) m i n « 
0.9). In the following, we therefore impose a more restric- 
tive jet-jet separation cut of AR(jj) > 1.0, which results in 
a sa 45% reduction of the background cross sections while 
reducing the signal only by about 7 — 8%. Note that in the 
large m t limit the AR(jj) ni i n distribution of the signal peaks 
at AR(jj) m i-a = and drastically differs in shape from that 
calculated using the exact loop matrix elements. If one were 
to calculate the signal cross section in the large m t limit, a 
AR(jj) cut would not result in a reduction of the background. 

Unfortunately, with two neutrinos present in the final state, 
\/l cannot be reconstructed. However, we anticipate that the 
invariant mass of all observed final state leptons and jets, m V i S , 
will retain most of the expected behavior of the different pro- 
duction processes. Figure || clearly demonstrates that this is 
the case: the signal peaks at lower values of m V i S than ei- 
ther background, especially for lower Higgs boson masses. 
However, the WWWjj background has a significant con- 
tribution from WH(-^ W + W~)jj production, resulting in 
a m V i S distribution which is similar in shape to that of the 
HH signal. Whereas the signal is concentrated in the region 
m V i S < 500 GeV, the background processes have a significant 
tail extending to m V i s = 1 TeV This makes it possible to nor- 
malize the background using data from the m v i s > 500 GeV 
region. Using exact loop matrix elements, the signal displays 
a pronounced peak which gradually moves to higher values of 
m V i S with increasing Higgs boson mass. In contrast, in the 
large m t limit, the m V i S distribution of the signal is extremely 
broad. 







pp^ ± f ± +4j 


- 




~____ 300 fb 1 


95% CL limits 




- 




600 fb" 1 




- 
- 


- 


3000 fb -1 




- 


- 600 fb -1 




SM 














300 fb" 1 







140 



160 



180 



(GeV) 



200 



95% 



CL 



for 



Figure 4. Limits achievable at 

AXhhh = Xhhh - 1 (Xhhh = X/Xsm) mpp — > t^l'^ + 4j at 
the LHC. Bounds are shown for integrated luminosities of 300 fb -1 
(solid lines), 600 fb -1 (dashed lines) and 3000 ftT 1 (dotted lines). 
The allowed region is between the two lines of equal texture. The 
Higgs boson self-coupling vanishes for AXhhh = — 1. 

The Feynman diagrams contributing to gg — > HH in the 
SM consist of fermion triangle and box diagrams @]. Non- 
standard Higgs boson self-couplings only affect the triangle 
diagrams with a Higgs boson exchanged in the s-channel. 
They only contribute to the J = partial wave, and thus 
impact the m v i s distribution mostly at small values. This is 
illustrated in Fig. |] for Mg = 180 GeV and two values of 
^hhh = A/AgM- Since box and triangle diagrams interfere 
destructively, the gg — > HH cross section for 1 < Xhhh < 
2.7 is smaller than in the SM. The absence of a Higgs boson 
self-coupling (Xhhh — 0) results in a Higgs pair production 
cross section which is about a factor 3 larger than the SM re- 
sult. 

The shape change of the m V i S distribution induced by non- 
standard values of Xhhh can be used to derive quantitative 
sensitivity bounds on the Higgs boson self-coupling. Figure^ 
shows the 95% confidence level (CL) limits for AXhhh = 
Xhhh — 1 which are obtained from a x 2 test of the m V i S 
distribution. The allowed region is between the two lines of 
equal texture. In deriving the bounds displayed, we combine 
channels with electrons and muons in the final state, conserva- 
tively assuming a common lepton identification efficiency of 
e = 0.85 for each lepton. In order to approximately take into 
account the (small) contributions to the background from titi, 
W Z + Aj and tij production which we have ignored in our 
analysis, we scale the background differential cross section by 
a factor 1.1. As mentioned before, our calculation of the back- 
ground processes does not include QCD corrections which are 
expected to modify the relevant cross sections by 20 — 40%. In 
order to derive realistic limits, we therefore allow for a normal- 
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ization uncertainty of 30% of the SM cross section. Since the 
background cross section can be directly determined from the 
high m V is region with a statistical precision of 15% or better 
for the assumed integrated luminosities, the bounds we derive 
are conservative. 

We derive sensitivity limits for integrated luminosities of 
300 fb -1 , 600 fb" 1 and 3000 fb"\ and Higgs boson masses 
in the range 150 GeV < M H < 200 GeV. Outside this range, 
the number of signal events is too small to yield meaning- 
ful bounds. For Mh < 150 GeV, this is due to the small 
H -> W*W branching ratio. For M H > 200 GeV, the 
gg — > HH cross section is too small. An integrated lu- 
minosity of 300 fb -1 (600 fb -1 ) corresponds to 3 years of 
running at the LHC design luminosity with one (two) de- 
tectors. The larger value of 3000 ft -1 can be achieved in 
about 3 years of running if the planned luminosity upgrade to 
C = 10 35 cm~ 2 s _1 [pj|] is realized. Figure Q demonstrates 
that, for 300 ft -1 , a vanishing Higgs boson self-coupling 
(AXhhh = —1) is excluded at the 95% CL or better, and 
that A can be determined with a precision of up to 60%. Dou- 
bling the integrated luminosity to 600 fb -1 improves the sen- 
sitivity by 10 - 25%. For 300 fb" 1 and 600 fb" 1 , the bounds 
for positive values of AXhhh are significantly weaker than 
those for AXhhh < 0, due to the limited number of signal 
events in this region of parameter space. For 3000 fb -1 , the 
Higgs boson self-coupling can be determined with an accuracy 
of 20 - 30% for 160 GeV < M H < 180 GeV. 

In summary, inclusive pair production of Higgs bosons at 
the LHC, with decays to a same-sign lepton pair and four 
jets via four W bosons, will make it possible to perform a 
first, albeit not very precise, measurement of the Higgs boson 
self-coupling A. The non- vanishing of A can be established 
at 95% CL or better for 150 GeV < M H < 200 GeV. The 
bounds on A derived here should be viewed as approximate. 
They can probably be strengthened by including other final 
states such as 3£ + 2j, or by using more powerful statistical 
tools than the simple \ 2 test we performed. More details of 
our analysis, along with inclusion of additional channels, will 
be presented elsewhere jl8[]. 
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